Casimir force in the presence of a medium 



Fardin Kheirandisffl 

Department of Physics, University of Isfahan, Isfahan 81746, Iran and 
Quantum Optics Research Group, University of Isfahan, Hezar Jarib Ave., Isfahan, Iran. 

Morteza Soltan{0 and Jalal Sarabadanjf] 

Department of Physics, University of Isfahan, Isfahan 81746, Iran 

In this article we investigate the Casimir effect in the presence of a medium by quantizing the 
Electromagnetic (EM) field in the presence of a magnetodielectric medium by using the path in- 
tegral formalism. For a given medium with definite electric and magnetic susceptibilities, explicit 
expressions for the Casimir force are obtained which are in agree with the original Casimir force 
between two conducting parallel plates immersed in the quantum electromagnetic vacuum. 

PACS numbers: 12.20.Ds, 03.70.+k, 42.50.Nn 



I. INTRODUCTION 

One of the most remarkable and fundamentally im- 
portant result of the field quantization is the Casimir 
effect which is a force arising from the change of the 
zero point energy caused by imposing the boundary con- 
ditions (BC) pp. This force is the macroscopic aspect 
of the quantum electrodynamics that provides a direct 
line between quantum field theory and the macroscopic 
world. The original calculation of the Casimir force be- 
tween two perfectly conducting parallel plates immersed 
in the quantum electromagnetic vacuum is based on the 
definition of the Casimir energy in the presence and the 
absence of boundary surfaces pQ that leads to an attrac- 
tive observable force 
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between the plates, where h is the Planck constant, c 
is the speed of light and H is the distance between the 
plates. In other way one can consider this effect by eval- 
uating the radiation pressure on macroscopic objects [2]. 
As the magnitude of the Casimir force is substantial at 
H < 100 nm this effect is relevant in nano-technology 
[6] and should be take into account to design and actuate 
microelectromechanical (MEMS) systems. Moreover the 
possibility of transducing the energy from the vacuum is 
investigated by means of MEMS [SHTU], 

Many attempts have been focused on observing the 
Casimir force and performing high-precision measure- 
ments during last few years [TT| - H^] . All these experi- 
ments are in agree with the prediction of the Casimir [1 
within a few percents. These deviation from the ideal 
force may due to temperature, roughness of surfaces and 
finite dielectric constants that have been covered in a very 
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recent book entitled by Advanced in the Casimir Effect 

m 

In 90th decade, Golestanian and Kardar developed 
a path integral approach to investigate the dynamic 
Casimir effect in the system of two corrugated conduct- 
ing plates surrounded by the quantum vacuum |181 119] . 
Emig and his colleagues also used the path integral for- 
malism to obtain normal and lateral Casimir force be- 
tween two sinusoidal corrugated perfect conductor sur- 
faces [30] [21] • Later on, the exact mechanical response 
of the quantum vacuum to the dynamic deformations of 
a cavity and the rate of dissipation have been calculated 
by using path integral scheme [52] . This motivated us to 
investigate the Casimir effect in the presence of a mag- 
netodielectric medium by quantizing the electromagnetic 
(EM) field using path integral formalism. Our system 
contains of a magnetodielectric medium with permitivity 
e and permeability /i, enclosed by two semi-infinite ideal 
metals (el — > oo for the ideal metal in left-hand side and 
£r — » oo for right-hand side one) as depicted in Fig. ([I]). 
We model the magnetodielectric medium by a continuum 
of harmonic oscillators (Hopfield Model) [2"5H2"5] . 

The outline of this paper is as follows: In Sec. [n] in or- 
der to introduce the scheme, we first quantize the simple 
case of a scaler Klein-Gordon field in the presence of a 




FIG. 1: This picture illustrates the schematic figure of the sys- 
tem under consideration. A magnetodielectric medium (e, fi) 
is enclosed between two perfect parallel conductors (er and 
el — > oo). The distance between conductors is H, and the z 
direction is perpendicular to the surfaces of the media. 
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medium. Sec. Mil is devoted to obtain the Casimir force 
for the scalar filed in the presence of a medium for differ- 
ent kinds of boundary conditions. In Sec. |IV| we develop 
our formalism to the case of EM field in the presence of 
a magnetodielectric medium. Sec. [V] gives the Casimir 
force for the case of EM field. Finally, the conclusions 
and outlooks are in the Sec. I VII 



II. 



FIELD QUANTIZATION USING PATH 
INTEGRALS 



To illustrate the method and also for later convenience, 
before considering the EM field in the presence of a 
medium, we consider the simplest case i.e. a scaler mass- 
less field. In the next section we will show that the Klein- 
Gordon field can be corresponded to each polarization of 
the EM field. Therefore, let us consider the following 
Lagrangian for the total system 



C = C S 



where 



C 



n/« - 2 ^fd^f, 



(2) 



(3) 



is the Lagrangian density of the massless Klein— Gordon 
filed. The medium is modeled by a continuum of har- 
monic oscillators as 126 
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d^i-pY^-p^Y*), (4) 



where Y^ is an oscillator's field, p is the density of mat- 
ter field and the interaction between the system and its 
medium is defined by 



where 



P 



C-int — fP, 



du)i/(ui)Y u 



(5) 



(6) 



In the next section we will show that the quantity P is in 
fact the polarization field corresponding to the medium 
and the interaction ^ will become the electric-dipole 
interaction. 

Generally a generating function is defined by [S7J 

Z[J] = J 2?[V>]exp ji j d n+1 x[£(tp(x)) + J(x)ip(x)] 

(7) 

where ip is the scalar field and the different correlation 
functions can be found by taking the repeated functional 
derivatives with respect to the source field J(x). The 
above partition function is Gaussian since the integrand 
has quadratic form with respect to the fields. To obtain 



the generating function for the interacting fields, we first 
calculate the generating function for the free fields 

z [j v ,j u ] = I v[<p\Y\v\r u ] x 

xexpjz J d n+1 x[C sys + £ mat + J v ip + J duiJ^Y^] 



(8) 



Using the n-dimensional version of Gauss's theorem we 
find 

J d n+1 x d^tp d^ip = - J d n+1 x ip □ if, (9) 

where □ is the d'Alemberian in (n+l)-dimensional space- 
time and the integration by part 



cT^xYuYu 
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d n+1 rY —Y 
at 1 



(10) 



the free generating function Q can be written as 
Zo [J*, J u ] - J Hv\Y u ]D[<p] exp | - % - J d n+ h 



<p(x)D(p(x)+ j duY u {x){— + pu> 2 )Y u (x) 



+ Jtp(x)ip(x) + / du)J u (x)Y u (x) 



(11) 



The integral in the equation (111 can be easily calculated 
from the field version of the quadratic integrals and the 
result is 



Zo [Jip, Jlj] 



exp 



d n+1 x I <r +L x 



71+1 



J v (x)G (x - x')J v (x') 



+ dcjJ 0J (x)G u (x — x')J u (x') 



(12) 



where Gu>(x — x') and Gq(x — x') are the propagators for 
free fields and satisfy the following equations 



UG (x-x') =S(x-x'), 



(13) 



{p— +pu J 2 }G LU (x-x') = S(x-x'). (14) 

We employ the Fourier transformation to solve the equa- 
tions (13) and (14). The solutions are 



Go(x-x') 



and 



1 



(2 7 r) n + 1 



d"kduj 



^k- (x— x )— iaj(t — t') 
UJ 2 -k 2 ' 



(15) 



1 f e -iu'{t-t') 

^-^=Ni dw W^- x ')' (16) 
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Here the space component of the point x £ lZ n+1 is indi- 
cated by the bold face x £ TZ n and the time component 
by t or Xo £ 1Z. 

For further use wc define 



Jp( z ) = / duv(u)J u (z), 



(17) 



the generating function of the interacting fields can be 
written in terms of the free generating function as [27] 



n ' - t/ 



n=0 



ipi 
1 n 



SJJz) dz SJ v (z) 



(18) 



where Z[Q] is the partition function of the free space. 
Thus the Green's function of Klein-Gordon filed can be 
obtained via 



G vv (x -y) = i 



5 2 Z[J ipi J P ] , 
8JJx)5JJyy J ^=°' 



(19) 



Combining the generating function (18), the Green's 
function ( 19 ) and the definition of C int (|5| yield the fol- 
lowing series for the Green's function 



G vtp (x - x') = G (x - x') 



O 2 



+ du> dx 1 dx 2 G (x-x 1 )i' (ui)—-^G ul (x 1 -x 2 )G (x 2 -x') 



dt 2 



i:) 2 



+l<kdldu)'ldxidx2ldx3dx4,Go(x—Xi)v (lo)--^G U j(xi—x 2 )x 



xG (x 2 -x 3 )v 2 (u')-^G u ,(x 3 -x$Go(x i -x')+- ■ ■ . (20) 

It is appropriate to use a general Green's function in (n + 
l)-dimensional Fourier space that is 



G(k, u) = J e lk - x ~ luJt G{x) dt d n x. (21) 

Therefore the Green's functions of the free fields ip and 
in the Fourier space might be given by 



G (k,w) 



1 



k 2 



(22) 



and 



G w (k,u/) 



1 fdw"dtd 3 x 



-tuj" t 



LU 2 — Uj" 2 — t0 + 
1 



p UJ 2 — uj' 2 



^0^ 



6(x)e 4((yt - k ' x) 

(23) 



respectively. Since we are interested in retarded Green's 
functions we have added —i0 + to the denominator of the 



equation (23). Since the reservoir field is assumed to 



be homogeneous, the Green's function of the reservoir 
does not depend on k in the above equation. Using the 



equations (22) and (23), G vv (x — x') can be written in 



the Fourier space as 
G IPilp (k,cj) 

= G (k, + J2i / dJuj 2 v 2 {J)G^ (w)Go(k, w)H 
Go(k,w) 



1 - fdu'u> 2 v 2 (u]')G a ,(Lj)G (k,u)) 
1 

V 2 -^ 2 - 1 - fdio' ,f ( ^ ) "n + ' 



(24) 



This Green's function can also be obtained directly from 
the Heisenberg equations of motion. By direct substitu- 
tion we can show that the Green's function G (fi(p (x — y) 
satisfies the equation 



□G w (x-x',t-t') 

_d_ f* ,.„ ,. d 
dt. 



dt"x(t-t") — G w (x-x',t"-t') = 6(x-x',t-t'), 

(25) 



which is the motion equation of dissipation field with sus- 
ceptibility of the medium x(i), with the following Fourier 
transform 



1 



y(u>) = — dui' — 77, 
Ay ' pj uj' 2 - uj 2 + 10+ 



(26) 



From the equations ( 24 ) and ( 26 ) it is clear that the mod- 
ified Green's function G (/3jV (k,w) can be obtained from 



the free field Green's function Go (a; — x') in Eq. ( 15 1 sim- 
ply by replacing uj 2 with e(lu)uj 2 , where e(w) := 1 +x(uj). 
It can be easily shown that this susceptibility satisfies the 
Kramers-Kronig relations as expected. 

By the same technique we can obtain correlation be- 
tween the polarization field and the Klein-Gordon field. 
To this end we define G v .p as 



G v ,p — 



S 2 Z[J V ,J P ] 
5J v 5Jp 



(27) 



that can be obtained via direct calculation similar to the 
procedure that ended to G Vilp (lc,u>) as 

G v .p{x - x ) = / dxiG (x - xi)--—G 0P (xi - x ) 



dz 



Q2 

+ I dx 1 dx 2 G p{x - a;i)^jG (xi - x 2 )G p(x 2 - x') 



+ j dx ± dx 2 J dx 3 dx4,G p{x ~ xi)-^G (xi ~ x 2 ) x 
d 2 

xG 0P (x 2 - X3)-^G {x a -x i )Gop(x 4 -x ) H ,(28) 



where 



Gop(x — x 1 ) = I dwv 2 (uj)G u {x — x 1 ). 



(29) 



4 



If we again write G lfi> p(x — x') in the Fourier space we 
find 

G Vt p(k, w) = iojGo(k,u})Gop(k,uj) 

bj 2 y 2 {w')G M i (w)Go(k, w)]" } 
G OP (k, w )Go(k, w ) 



1 -/rfw'w 2 ^ 2 (w')G^(w)Go(k,u;) 
G p(k,w) 



k 2 -cj 2 - i fdw'-^ 

p J UJ' 

Comparing Eqs. ( 16 ) and ((29 

J uj ■ 



v 2 (oj')oj 2 



(30) 



rewritten as 



+i0+ 

yields Gop(k, uj) = 
+i0+ — x(w), consequently G Vi p(k,w) can be 



G V: p(k, uj) — iujx{uj)G v ^{\t, u}), (31) 

The other important correlation function is Gp.p{x — x') 
which is defined via generating function Z[J V , Jp] as 



Gp, p (z-0 = 



S 2 Z[J V ,J P ] 



SJpSJp 

By straightforward calculations one 
Gpp(k, uj) as 



(32) 

can obtain 



Gp,p(k, uj) = 



v 2 {uj) 



+ uj 2 x 2 (uj)G lpip (k,uj). 



(33) 



The imag inar y part of the response function can be read 

I I 2 / \ 

from Eq. (26) as \ ' — Jm^(w). Here to illustrate the 



validity ot our results, i.e. Eq. (31 33 1, we compare them 



with the results of the other methods of field quantiza- 
tion. In other conventional methods of phcnomcnological 
field quantization |28) , the fields can be divided into posi- 
tive (+) and negative (— ) frequencies parts which satisfy 
the constitutive relation 

P ± (k,uj) =±iuj X (u J )tp ± (k,u J ) + P N (k ± ,to), (34) 

after quantization of the fields, where by "we mean oper- 
ator and the operator with positive frequency is the Her- 
mitian conjugate of the negative one. P N is the noise part 
of the polarization field that according to the fluctuation- 
dissipation theorem [25] 128] satisfies 



[P#(x,w),P£(x,w)] =7r/mxM«5(x-x'), (35) 

which [• • • , • • • ] denotes the commutator of two opera- 
tors. If we use the Eqs. (34) and (35) to obtain the 



Green's functions, we achieve the same results as the 
Eqs. ( 3l][33 ). This shows the validity of our path inte- 
gral quantization. 

According to the definition of the Green's functions, 
after integrating over gaussian fields one can read the 
generating function ( 18 1 as 



Z[J v ,Jp] 



- exp^i 
-JipGp^ 



dx J dx 1 [J^G^ 



x ) J \p 



x — X 



')J P + J P Gp.p(x - x')J P ] 



(36) 



III. CALCULATING THE CASIMIR FORCE 

A. General formalism 

In this section we briefly review the path integral tech- 
nique to calculate the Casimir force. Let us consider two 
conducting plates faced each other at the distance H and 
embedded in an arbitrary medium. The field (p satisfies 
the Dirichlet 



or Neumann 



<p(X a ) = 0, 



d n ip(X a ) = 0, 



(37) 



(38) 



boundary conditions on surface, where X a , (a = 1,2) 
is an arbitrary point on the ath conducting plate. To 
obtain the partition function from the Lagrangian we use 
the Wick's rotation, (t —> it) and change the signature 
of the space-time from Minkowski to Euclidean. The 
Diriclet or Neumann boundary conditions can be taken 
into account using the auxiliary fields ip a (X a ) |22) 

5(<p{X a )) = f Vty a (X a )yf dX ^ x ^ x «\ (39) 



and 



S(dMX a )) = / V[^ a (X a )]e-^ dx ^ x ^ x -K 



(40) 

After Wick's rotation the Dirichlet and Neumann parti- 
tion functions can be cast into the form 



and 



r 2 

Z D = Z^ 1 / T> [(p] J^J T> [ip a (X a )] ) e s ° 

J a = l 

Z N = Z^ 1 [ V[<p}\{V^ a (X a )])e 

J a = l 



(41) 



(42) 



respectively, where Zq is the partition function of the free 
space, and 

S D [p] = { d {n+1) x 



{C(ip(x)) +tp(xj$2 d^XS(X-X a )^ a (x)}.(43) 



and 



S N [cp] = J d {n+1) x 
x{£^(x))+^(x)Y^ I d^XS(X - X a )d n ^ a {x)}.{U) 
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Using the same procedure of Ref. ( |22j ) the parttion func- 
tions for the Dirichlet and Neumann BC can be read 



Z D = 



1 



and 

where 

T D (x,y,H) 

and 

T N (x,y,H) = 



\J det Tn(x, y, H) ' 



g(x-y,Q) G(x-y,H) 
G(x-y,H) G(x-y,0) 



(45) 



(46) 



(47) 



-d 2 z g(x-y,0) ~d 2 z G(x-y,H) 
-d 2 z g(x-y,H) -d*g(x- y ,o) 



(48) 

where g is the Green's function of the fields after wick 
rotation. We define the effective action as 



-ilnZ(H), 



(49) 



where In Z(H) can be either for the Dirichlet or Neumann 
BC, in order to calculate the Casimir force by applying 
derivative with respect to the distance between the plates 



F 



dSeff(H) 

dH ' 



(50) 



It is easy to show that the contribution of the Dirichlet 
is the same as that of the Neumann BC to the Casimir 
energy and hence the Casimir force in the presence of a 
isotropic and homogenous medium like [21 j . So that in 
the next sections we treat only the Dirichlet BC. 



B. Casimir force for different boundary conditions 

In this section we would like to obtain the Casimir 
force in the presence of an absorptive medium. Before 
we obtain the Casimir force for interacting fields, we in- 
vestigate the possibility of the existence of the Casimir 
force due to the matter field alone. Using the Lagrangian 
Q and expression for the effective action ( 49 ) we find the 
Trj tensor as 



T u (x,y,H) = 



g u (x 
g^{x - 



-y,o) g^ix- 
y,H) <L(x 



-v,o) 



(51) 



But since Gu(x — y, H) = for this situation, the nonin- 
teracting matter field alone, does not lead to any modi- 
fied Casimir force. This result is clear since we model the 
matter field by the Hopfield model ****** [HOPFIELD's 
REFERENCE]******. This model is based on an in- 
dependent set of harmonic oscillators and imposing any 
condition on one of these oscillators does not affect the 
others, and hence we do not expect any Casimir effect. 



For a dissipative field (p (251 we may consider three 



different boundary conditions, 

i) Imposing the boundary condition on the Klein- 
Gordon field: for this case the T tensor can be read 



r vv (x,y,H) 



G<p,<p(x — y,Q) g VlV {x -y, H) 
g v ,v(x-y,H) G v ,<p(x - y,Q) 



(52) 



Since is diagonal in the Fourier space, to obtain 
the Casimir force we proceed in this space. The Fourier 
transformation of G vv (x — y, H) is 



g vv {p,q,H) 



J dxdye lp x+tq y G vv (x-y,H) 



e -n(p )\po\h 

2n(Po)bo 



(53) 



where p = (po, p), p is a vector parallel to the conductor, 
Po the temporal component of the p, n(po) = y/ s{po) and 
e(po) = e(iui). Thus for the case i the Casimir force is 



Fi = - 



d 3 P , E(p) 



(54) 



where E{p) — [n 2 (po)Po + p 2 ] 1 ^ 2 - In the absence of the 
medium between the conductors, n(po) — 1, we recover 
the original Casimir force between two plates immersed 
in the quantum vacuum of a scalar field 



F t = - 



d 3 p 



( 27r )3 e 2\ P \H _ j A80H 4 ' 



(55) 



and for a non absorptive medium with the susceptibility 
x(t) = XoS(t), we find the modified Casimir force as 



Fi — —Fy ac . 
n 



(56) 



The above relation is fully in agree with the result of 
the Lifshitz theory of fluctuation-induced force bewteen 
media 30J , [31 . The equation (54) is interesting since it 
is the reminiscent of the Bose-Einstein distribution. In 
fact, E(p) can be interpreted as the force density due to 
the bosons in the state p. 

ii) Imposing the boundary condition on the polariza- 
tion field: in this case Tpp tensor is 



T PP (x,y,H) = 



gp,p(x-y,0) gp,p(x-y,H) 
gp,p(x-y,H) gp,p(x-y,0) 



Here the Casimir force is 
d 3 p 



Fa = - 



[X 2 (P0] 



E(p) 



ae 2E(p)H _ X J 



(57) 



(58) 



where a = E(p)Imx(po) + X 2 (Po)- Although the noise 
operators do not have any spatial correlation but their 
presence on the surface can affect the Casimir force and 
decrease it due to polarization. 



G 



iii) Imposing the boundary condition on the both of 
polarization and Klein-Gordon fields: we can easily show 
that T VVt pp(x, y, H) is a 8 x 8 tensor with the form of 



T vlp (p,q,H) 



9oX(<Zo)r w O,g,H) 



gox(?o)r w O, QiH) %X (2o)IW(p, q,H) 



-Imx(qa)I 
(59) 



where I is the 2x2 unit matrix multiplied by (27r) 3 <5(p + 
q). It can be easily shown in this case that the Casimir 
force will be the same as that of imposing the bound- 
ary condition only on the Klein-Gordon field i.e. case i. 



We can interpret this situation with the aid of Eq. ( 34 1 



According to this relation if the BC is imposed only on 
the Klein-Gordon field then is not zero, and if the 
BCs are imposed on the both of polarization and Klein- 
Gordon fields then the BC will be imposed on PS au- 
tomatically which according to the begining of the Sec. 
|IIIB| does not lead to aany Casimir effect. 

The physical interpretation of the first and third con- 
dition is obvious. These conditions arise when we want 
to calculate the Casimir force between two perfect con- 
ductors that enclose a medium. But the second BC can 
happen when we want to consider the system of contain- 
ing two dielectric slabs. This kind of boundary condition 
leads to the Casimir force between two dielectric slabs 
which is under consideration. Here we only consider the 
boundary conditions in cases i and iii, and we shall treat 
the case ii in elsewhere 1291. 



IV. ELECTROMAGNETIC FIELD 
QUANTIZATION IN THE PRESENCE OF A 
MAGNETODIELECTRIC MEDIUM 

In this section we develop the formalism to EM field 
in the presence of a magnetodielectric medium [35] . The 
Lagrangian of EM field in the presence of a medium can 
be written as 



£ = £ 



EM 



Co 



Ci 



where Lem is 



C 



EM 



(60) 



(61) 



where E and B are the electric and magnetic fields. They 
can be written in terms of scalar and vector potentials U 
and A respectively as E = A — VU and B = V x A. In 
this work we use the Coulomb gauge V • A = 0, i.e. A 
is a transverse field. Li ma t and £2mat reffer to the po- 
larization and magnetization of the medium respectively 
and can be written as 



i 



2 xL), 



(62) 



where Xj w is an oscillator's vector field. The interaction 
part of the Lagrangian is 

C mt = A P VC/.P + V x A ■ M, (63) 

where P and M are polarization and magnetization of 
the medium defined by 



/•OO 

P = / i^(w)Xi u 
Jo 

f- DO 

M = / ^ 2 (u;)X 2li 
Jo 



(64) 



The Euler-Lagrange equations lead us to the fact that U 
is not an independent dynamical variable. To obtain U 
in terms of the other independent dynamical variables, it 
is appropriate to write the fields in Fourier space where 
the longitudinal and transverse parts of a field can be 
separated using the unit vectors ea(k) = k and e^(k) 
(A = 1,2) respectively. ei(k) and e 2 (k) are perpendic- 
ular to each other and k. Using these unit vectors the 
scalar potential can be written in terms of the longitudi- 
nal part of the matter field as (In what follows we indicate 
the fields in Fourier space by a ~ over them.) 



U 



£o|k| 



(65) 



By applying this recent relation to the Eq. (601, we can 
easily show that the longitudinal part of the electromag- 
netic field only changes the longitudinal component of 

£lmat aS 



f\\ _ 
l -lmat ~ 



./2v\\2 



(66) 



where u>' = y/w 2 + (J 2 and lo 2 = v . It is worth noting 
that only the transverse parts of Lagrangians have con- 
tribution to the Casimir effect and the longitudinal part 
does not lead to any Casimir force. Therefore we just 
consider the transverse part of the Lagrangian which in 
the Fourier space is 



L 



where 



£i m = e (A 2 -c 2 B 2 ), 



(67) 



(68) 



AU* = / dLu(p±t 2 - P uj 2 ±t 2 ), (69) 



and C-i int an d ^2 int re f er to the interaction of the polar- 
ization and magnetization fields with the electromagnetic 
field respectively which are 



(70) 
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and 



._>_,:„, , dwi/ 2 (w)k x A ■ 



9 

^ / dw^ 2 (w)|k|i A X 2 > 



£AV, (71) 



where ew is the antisymmetric tensor. 

To obtain the Casimir force between two plates with 
axial symmetry due to the transverse part of the La- 
grangian, the filed modes can be divided into TM and TE 
modes [501 HI] • It can t> e shown that for our case as we 
consider the electromagnetic field in the presence of a ho- 
mogenous, isotropic and flat magnctodiclcctric medium, 
similar to [20, 21j again the field can be divided into TM 
and TE modes which refer to the Dirichlet and Neumann 
BC respectively. As for the situation under study here 
for TM and TE modes the Casimir force is the same, 
therefore here we consider only the TM mode. Accord- 
ing to the above discussion the full Lagrangian can be 
rewritten as 



C = C S 



where 



C-sys = \ d^(pdfj,(p, 



(72) 



(73) 



is the Lagrangian density of a massless Klein— Gordon 
filed and 



£mat — 



2 

E 



i • i 

du;(-pXt--pu J 2 Xl). 



The interaction term is defined by 



(74) 



(75) 



where P = J and M = J duJV2(oj)X2 U - It 

can be seen that the Lagrangian ( |72| ) is similar to the 
Lagrangian (pi). In fact these Lagrangians are the same 
if we take v^ytS) = 0. This is the reason why we called 
P as a polarization field and the results obtained in the 
Sec. |n] can be used for a polarizable medium. 

By mixing (|18[) and ( 75 ) , and use the same procedure 



in the Sec. |II[ after some manipulations we obtain the 
Green's function as 



G Vtlp (k,(jj) = 



1 



k2(l-XmH)-W a (l + X e H)' 



(76) 



which is the Green's function of EM field in the pres- 
ence of magnetodielectric medium with the electric and 
magnetic susceptibilities 



f + OC 



fiV) 



and 



Xm(w) = 



+ 00 



dJ- 



4W) 

-uj' + z0+ ' 



(77) 



(78) 



respectively. These susceptibilities satisfy the Kramers 
-Kronig relations as expected. The other Green's func- 
tions or correlation functions are 



(79) 



GL )A f(k,a;) = t\\t\wXm(u)G<p tl p(b,u)), (80) 



Gp,p(k,w) 



+ ^XeMG w (k, w ), (81) 



+ |k| 2 x ^( w )G w (k, w ). (82) 



The first terms in the right hand side of Eqs. ( 81 1 and ( 82 ) 



are related to the noise operators of the system which sat- 
isfy the fluctuation-dissipation theorem, like (35). Con- 
sequently the generating function is 



Z[J v ,Jp] — exp^i J dx J dx' 



J tp\pj)G tn^(n\X X )o^^(x) 



+J v {x)Gp. v (x - x')Jp(x) + J v {x)Gp tV {x - x')Jp{x) 



(83) 



This relation can be used to obtain the Casimir force in 
the next section. 



V. CALCULATING THE CASIMIR FORCE FOR 
EM FIELD IN THE PRESENCE OF A 
MAGNETODIELECTRIC MEDIUM 

Here again we consider three different boundary con- 
ditions similar to the Sec. Mil varies kinds of the fields. 

i) Imposing the boundary condition on EM field: if 
we impose the Dirichlet BC on EM field the partition 
function becomes 



Zn = 



1 



y/det T tptp (x,y,Hy 



(84) 



where 

T vv (x,y,hH) 



Q<p,<p(x - y,0) G v , v (x - y,H) 
G v , v (x - y,H) G v , v (x-y,0) 



. (85) 

and Gtp.tp {x — y,h) is the Green's function ( 76 1 with imag- 
inary time. 

To calculate Zp), we invoke the Fourier space where 
r w (i, y, hH) is diagonal and its elements have the form 

Q VtV (p,q,H) 

e -V" 2 (Po)Po+P 2 H 

= KPo h r^r-, 2 (2n) 3 S(p + q), (86) 



where here n(po) — \J m(po)e(po) and niuf) 
Finally the Casimir force is obtained as 



F 



d 3 p E(p) 



(2n) 3l e 2E W h - 1 J ' 



(87) 



where E{p) — \/n 2 (po)p 2 ] + p 2 . This equation is like the 
equation ( |54[ ) the only difference is in the definition of 
n(po). If we impose the Neumann BC on the EM field 
we will achieve the same result as that of the Dirichlet 
BC. 

ii) We can impose the BC on polarization and magne- 
tization fields. As we mentioned after the equation (59 1, 



this situation does not appear in our problem since we 
consider a magnetodielectric medium surrounded by two 
perfect conductors which lead to the BC on EM field. 
But for two magnetodielectric slabs this kind of BC may 



appear which may be dealt with elsewhere [25] , 

iii) We can impose the boundary condition on both 
matter and EM fields. In this case if we redo the calcu- 
lations, we conclude that this situation leads to the same 
result as case i. This shows that for a magnetodielectric 
medium the conditions i and iii lead to the same result 
for the Casimir force. 



VI. CONCLUSION 

In this article we quantized the electromagnetic field 
in the presence of a magnetodielectric medium in the 
frame work of path integrals. For a medium with a given 
susceptibility, the modified Casimir force is obtained for 
different boundary conditions. The present approach can 
be generalized to the case of rough perfect conductors in 
the presence of a general medium straightforwardly [29j . 
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